Ferromagnetic frustrated spin systems on the square lattice: a Schwinger boson study 
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We study a ferromagnetic Heisenberg spin system on the square lattice, with nearest neighbors 
interaction Ji frustrated by second J2 and third J3 neighbors antiferromagnetic interactions, using 
a mean field theory for the Schwinger boson representation of spins. For J3 = we find that the 
boundary between the ferromagnetic and the collinear classical phases shifts to smaller values of J2 
when quantum fluctuations are included. Along the line J2/IJ1I = 1 the boundaries between the 
collinear and incommensurate regions are strongly shifted to larger values with respect to the classical 
case. We do not find clear evidence for spin gapped phases within the present approximation. 

PACS numbers: 75.10.Jm; 75.30.Kz; 75.30.Ds. 



I. INTRODUCTION 

Most current works on frustrated magnetic systems 
generally deal with competing antiferromagnetic in- 
teractions. Recently, some frustrated systems have 
been discovered where the basic interaction is ferro- 
magnetic. In particular some vanadate and cuprate 
crystals as {CuCl)LaNb207^ Pb2VO{P04)2r'- 
SrZnVO{P04)2r- BaCdVO{POi)2^ and 
PhZnVO{P04)2^ can be described by a two di- 
mensional Heisenberg model of spin S ~ 1/2 with a 
ferromagnetic first neighbors interaction and antiferro- 
magnetic further neighbors interactions. Other possible 
relevant materials arc {CuBr)LaNb20-r^ which shows 
collinear order, and {CuBr)Sr2Nh^Oi[r^ which shows a 
plateau at M = 1/3 in the magnetization curve. As one 
dimensional counterparts, materials like LiCuVO^ir^^ 
and Li2ZrCuO^ can be modeled by ferromagnetic 
frustrated spin S = 1/2 Heisenberg chains. 

In the present work we consider such a two dimensional 
Heisenberg model on the square lattice (see Fig. [1} with 
ferromagnetic nearest neighbors interactions Ji < 0, frus- 
trated by next to nearest neighbors antiferromagnetic in- 
teractions J2 > and also third neighbors antiferromag- 
netic interactions J3 > 0, given by the Hamiltonian 

H = Ji Si.Sj + J2 Si ■ Sj 

+J3 J2 • ^1 (1-1) 

where Si is the spin S operator at site i and the range of 
the interacting neighbor sites i,j is indicated by brackets 
{i,j)r with r = 1,2,3. Closely related work has been 
done for antiferromagnetic JiJ^"— 

The classical S* — > 00 counterpart of these interactions 
is described by the scalar product of commuting vectors, 
where the lowest energy configuration is elementarily ob- 
tained. At any value of the exchange constants, it is 
described by a planar arrangement of vectors rotated by 



FIG. 1. (color online) Square lattice with Ji < ferromag- 
netic first neighbors couplings, J2 > antiferromagnetic sec- 
ond neighbors couplings and J3 > antiferromagnetic third 
neighbors couplings. 



relative angles "dx in x direction and i!)y in y direction, 
giving rise to the classical phase diagram in Fig. [5] com- 
posed of four different ordered phases fi2i2£ 

• F: a ferromagnetic phase, {■dx,'&y)={0,0). 

• CAF: a collinear antiferromagnetic phase show- 
ing antiferromagnetic order in one direction of the 
lattice and ferromagnetic order in the other one, 
{'&.,^y)= (0,^) or (^,0). 

• CH: a collinear helicoidal phase showing helicoidal 
order in one direction of the lattice and ferromag- 
netic order in the other one, {'dx,'&y)= (0, 9) or 
(g,0) with cos((7) = ^^ijp. 

• H: a helicoidal phase composed by helicoidal order 
in both directions of the lattice, {'dx,'&y)~{Q,Q) 
with cos(Q) = 2(j';+2j,) ■ 

Opposite to the classical limit, the quantum case S = 
1/2 has been recently analyzed using exact diagonaliza- 
tion (ED) techniques to explore the complete phase di- 
agram in Ref. [2l[, while more detailed features of the 
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FIG. 2. Classical phase diagram. F: ferromagnetic phase, 
CAF: coUinear antiferromagnetic phase, CH: coUinear heli- 
coidal phase, H: helicoidal phase. A vertical axis at J2 = | Ji 
is drawn for comparison with Fig. [31 



J3 = line were studied by ED in Ref. |22| - |24l | (and also 
by coupled cluster methods in Ref. H^l)- However, state 
of the art ED computations can reach system sizes only 
up to 40 sites. Moreover, discrepancies pointed out in 
Ref. [l^l with earlier results call for investigating these 
systems with complementary techniques. 

Motivated by the above discussion, we are interested 
in the effect of quantum fluctuations as introduced by 
large S methods. In a first step we have computed linear 
spin wave corrections to the classical order, as derived 
from the Holstein-Primakov bosonic representation for 
spins. While this method provides significant shifts for 
the classical antiferromagnetic Ji — J2 — J3 modelfi^ it 
gives very small corrections in the present ferromagnetic 
case (not reported here). However, it helped us to localize 
two particular areas in the phase diagram that seem to 
be clearly modified by the quantum fluctuations. These 
are: (i) the region of the classical phase transition be- 
tween the ferromagnetic (F) and coUinear antiferromag- 
netic (CAF) phases located at J3=0 and J2 « 0.5|Ji|, 
where ED studies show discrepancies;^ and (ii) the re- 
gion J2 ~ |Ji|, J3 > 0, where appcarcncc of a gapped 
phase was reported 1^ 

In the present work we analyze these two regions by 
means of the Schwinger boson (SB) representation for 
spins. This method does not start from any classical 
order (in contrast with linear spin wave theory) and, 
treated at mean field level (SBMFT), allows us to study 
quantum fluctuations in fairly large systems. Moreover, 
it has been tested to give quite accurate results even for 
S = 1/2, by comparison with ED (see e.g. Refs. (25l.[26|). 
Within SB theory, long range order is characterized by 
boson condensation. We classify the different ordered 
phases according to the condensate momentum and the 
evaluation of spin correlation functions, finding impor- 
tant shifts in the classical phase boundaries. A schematic 
phase diagram with our results is shown in Fig. [31 

The paper is organized as follows: in section [Hi we 
present the SBMFT methods, in section IIIII we analyze 
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FIG. 3. Corrections to the classical phase diagram computed 
from Schwinger boson mean field fluctuations. On the line 



J3 = we find the CAF phase for J2 > J2 



0.41 1 Jll 



for J2 > JI" = 0.58|Ji| self consistent solutions are obtained 
for systems up to 10* sites, while smaller system sizes are 
reached as J2 approaches J2~^^^ (empty marks). On the 
line J2 = |Ji| phase boundaries are strongly shifted to larger 
values with respect to the classical case. 



the model without J3 interactions, and in section ITVl we 
analyze the J2 = |Ji| line with J3 interactions. Section 
[Vlis devoted to the conclusions. 



II. SCHWINGER BOSON MEAN FIELD 
THEORY 

The Schwinger boson approach allows to incorporate 
quantum fluctuations while keeping the rotational invari- 
ance of the Heiscnbcrg model (see for instance Ref. ^27!]). 
In this method the spin operators are written in terms of 
two species of bosons and via the relations 



(2.1) 



In order to represent spin S properly, one must locally fix 
the bosonic occupation to 25-1-1 states by the constraints 
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(2.2) 



at each site i. 

The Heisenberg Hamiltonian is then a quartic form in 
bosons, but can be conveniently written as quadratic in 
bond operators, namely quadratic bosonic operators in- 
cluding one boson from each of the interacting bond sites. 
Such a factorization is not unique, and different schemes 
are adopted in case of ferromagnetio^ or antiferromag- 
netic frustrated^S interactions. A mixed schcm e^^'^° has 
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been shown better adapted to include both antiferromag- 
netic and ferromagnetic short range correlations. It deals 
with more mean field parameters, but provides quantita- 
tively better results and is our choice to study the present 
ferromagnetic frustrated system. 

Bond operators A and B are defined as 



(2.3) 



Notice that 
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where : O : means the bosonic normal order of an opera- 
tor O, relate non- vanishing A to antifcrromagnctic struc- 
tures and non-vanishing B to ferromagnetic structures. 
Moreover, expanding the squares yields representations 
for the SU{2) invariant terms Si ■ Sj. The Hamiltonian 
p.ip can then be written as 

H= E 'J^T.(-- ^h^^^i -Al.^^..') + Bx (2.5) 

'■=1,2,3 (jj), 

where the term 

Hx=J2 ^K&It^a + 4,M - 2S) (2.6) 

i 

forces the local constraints. A,; being the Lagrange mul- 
tipliers. 

At mean field level, we perform a Hartree-Fock decou- 
pling introducing a uniform Lagrange multiplier A and 
translationally invariant parameters a, /3 for the expec- 
tation values of each type of bond operator present in 
the Hamiltonian (as mentioned above, this decoupling 
is not unique). As it is known ;^^!^^" — the most severe 
approximation here is the violation of the local boson 
number constraint in Eq. (|2.2p . which is only respected 
on average. We are thus dealing with the Lagrange mul- 
tiplier A, six a's and six /3's as independent variational 
parameters, the latter set as expectation values of bond 
operators 



(2.7) 



and similar expresions relating the /3's to (Bij) expecta- 
tion values. For compact notation we write a^j = {Aij), 
Pi J = {Bi_j), using site indices to indicate the range of 
the bond {i,j)r {r ~ 1,2,3) as well as the possible ori- 
entations along the lattice described in Eq. (|2.7|) . The 



mean field Hamiltonian then reads 
Hmf 



= 1,2,3 (jjV 



Bh A,. 



a] 



■ E E (1/3^, 

r=l,2,3 {^,J)^ 



Hx.{2.8) 



After a Fourier transform one gets momentum modes 
with quadratic terms which are not particle number con- 
serving. These are diagonalized by a standard Bogoli- 
ubov transformation, depending on the variational pa- 
rameters and rendering decoupled modes with simple 
particle number conserving, positive, quadratic, Hamil- 
tonian 



Hmf = E 



-I- const, 



(2.9) 



where 




(2.10) 



contains the Bogoliubov bosonic operators, with disper- 
sion relation io{k), and const stands for non operator 
terms. 

Finally, we compute self consistently the mean field 
parameters by minimizing the ground state (Bogoliubov 
vacuum) energy with respect to A and equating aij and 
I3i_j with the ground state expectation values of the cor- 
responding operators. Such computation is done numer- 
ically on finite lattices of N sites with periodic boundary 
conditions, allowing the study of large system sizes, up 
to 10^ sites in the present work. 

We must stress that our procedure is not suited for the 
ferromagnetic phase, where parameters aij vanish and 
the Hamiltonian in Eq. (|2.8p is already particle conserv- 
ing: the Schwinger bosons vacuum simply violates the 
constraint in Eq. (|2.2p . even on average. We use then the 
exact energy of a fully polarized (ferromagnetic) state, 
Ep = 2NS'^{Ji + J2 + J3), for comparison with SBMFT 
energies, or extrapolations thereof, to determine the fer- 
romagnetic phase boundaries. 

Once the self consistent equations are solved, the tools 
above allow to compute any kind of observable on the 
ground state. In the present work we have set S = 1/2 
and studied four quantities: the dispersion relation, its 
gap, the modulated magnetization M,j (defined below) 
and the spin correlation function. 

When the dispersion relation shows a zero mode, Bose 
condensation indicates an ordered phase, in the sense 
that the spin structure factor shows a maximum at a 
pitch angle 9 = {9x,0y) commensurate with the finite 
lattice, related to the position of the zero mode of the dis- 
persion relation kmin by = 2kmin^ Notice that Bose 
condensation depends on boson density, related in the SB 
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approach to the spin S representation by the constraint 
in Eq. (|2.2p . As we study numerically the lowest density 
case, S = 1/2, such ordered phases will also be present 
for larger S. For illustration purpose, in Fig. 2] we show 
the dispersion relation at coupling values J2 = 0.75 |Ji| 
and J3 = (well inside the CAF classical phase). 
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FIG. 4. (color online) Dispersion relation for J2 = 0.75|Ji 
and J3=0 and for a system of size A'^ = 100 x 100 (in arbitrary 
scale, with darker zones indicating lower energy). The bosons 
condense at points kx=±n/2, ky=0, which correspond to a 
CAF phase (±7r,0). 

In the case of long range order, the vanishing of the 
gap in the thermodynamical limit is usually recovered 
after a finite size scaling analysis. Another issue arising 
at finite sizes is that related to commensurability. As 
the gap is obtained through the value of the minima of 
the dispersion relation on the reciprocal lattice, in the 
case that a thermodynamical minimum does not fit with 
the available momenta values at finite sizes, numerical 
difficulties may show up (see section ITVl . 

Regarding the other observables, the spin correlations 
{Si ■ Sj) for any pair of sites i,j are computed on the 
ground state from the bosonic representation in Eq. (|2.ip 
and the Bogoliubov transformation obtained numerically. 
Wc measure modulated magnetizations by considering 
the 6'-dependent susceptibility Af^(0), defined as the fol- 
lowing average over the lattice: 

^^"(^^) = jv(A2) ^^^-^^-^"""^''^'"''""^'^' ^^-^^^ 

where 6 = [Ox^Oy] and Ri is the position of spin Si 
with respect to some reference site. It amounts to an 
improvcdccl~ normalization of the spin structure factor, 
that fits better small systems and tends to moderate the 
weight of strong on-site terms. It is straightforwardly 
computed from the spin correlations. 

III. NEAREST NEIGHBORS FRUSTRATION J2 

In this section we analyze the case J3 ~ 0, that is a 
system with ferromagnetic first neighbors couplings and 



only second neighbors antiferromagnetic interactions. 

The classical phases on this line, shown in Fig. [21 are F 
and CAF, separated by a critical value J^"^"-""" = 0.5| Ji|. 
The quantum case was studied for S = 1/2: based on 
ED of the model and coupled cluster methods, Ricliter 
et. aB^ predict a simple shift of the critical coupling to 
lower J2 = 0.39| Ji|, while Shannon et al^^^ estimate by 
ED a CAF phase only for larger J2 > 0.6| Ji| and predict 
the presence of a quadrupolar (bond-nematic) phase in 
the critical area, 0.4 < J2/I Ji| ^ 0.6. 

Wc have studied with SBMFT systems of sizes ranging 
from = 4 X 8 (finding excellent agreement with ED, 
for ground state energies) up to A^=100xl00. 

As mentioned in Sect. [TTl our procedure does not pro- 
vide a self-consistent solution for the F phase. We first 
analyze the values of J2 above which SBMFT solutions 
are obtained. These values turn out to be sensitively 
dependent on the system size. Above J2 ~ 0.56|Ji| we 
reach solutions for all explored systems, up to 100 x 100 
sites; but approaching the F phase we get oscillatory be- 
haviour with the system trapped in metastable configu- 
rations, and the tractable sizes reduce as down as 20 x 20 
at J2 « 0.4|Ji|. The size dependence of the lowest cou- 
plings J2{N) tractable within SBMFT is roughly linear 
in 1/N, as shown in Fig. [51 suggesting an infinite size ex- 
trapolation to J2°° = 0.58 I Ji|. Thus wc estimate that, 
investing enough CPU time, one can treat systems of ar- 
bitrary size only when J2 > J2°°. 
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FIG. 5. Evolution of the lowest couplings J2{N) tractable 
within SBMFT along the line J3—O with the inverse of the 
system size. Temptative linear extrapolation realized for sizes 
TV x iV with iV = 52, 56, 60, 68, 80, 100. 

For J2 above J|° the observables computed from the 
SBMFT self-consistent solutions correspond to a CAF 
phase, showing staggered magnetization along one of the 
lattice axes. An example of the dispersion relation, for 
J2 = 0.75 I Ji I and J^—O in a large system of size N ~ 
100 X 100, is the one shown in Fig. [H The boson modes 
become gapless at momentum points k = (±7r/2,0), 
showing that the bosons do condense. The condensation 
momenta correspond to ordering angles 9 = (±7r, 0). The 
same pattern (alternatively with 6 — (0, ±7r)) is found for 
J2 > J2{N),N = 20,40,60,80,100. 

The dispersion relation gap goes clearly to zero for 
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A'' ^ oo, as shown in Fig. |6] (size scaling for J2 = | Ji| 
is shown in the inset). Correspondingly, the spin cor- 
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FIG. 8. (color online) Staggered magnetization as order pa- 
rameter for the CAP phase, at J3 = 0. 



FIG. 6. (color online) Energy gaps (minima of the dispersion 
relation) for systems of different size, in the CAP phase. Ex- 
trapolated value of the gap vanishes, confirming collinear an- 
tiferromagnetic order. Inset: gap extrapolation for J2 = |Ji|. 

relation function exhibits long range order: when 6 = 
(±7r,0), we observe antiferromagnetic correlations in the 
X direction and ferromagnetic correlations in the y di- 
rection. For example, in Fig. [7] we show the correlations 
for J2 = 0.75|Ji|. The corresponding modulated mag- 
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indications of the CAF phase. It turns out that the fi- 
nite size available solutions scale in this range to the same 
CAF phase as in the range J2°° < J2 < \Ji\- 

The energy per site obtained from SBMFT self- 
consistent solutions is very stable against system size and 
shows a neat linear dependence with J2, including the 
region 0.4 < J2/IJ1I < 0.58. (see Fig. I^. We conclude 
that lack of convergence of SBMFT self-consistent equa- 
tions at large system size is an artifact of our present 
approach, presumably due to the proximity of a ferro- 
magnetic phase, and that solutions obtained for small 
systems are good estimates of the CAF phase in the 
region for 0.4 5. J2 / Ui I < 0.58. Then, following the 
criteria in Ref. [23| we set the F-CAF phase boundary 
at the intersection point between the extrapolated CAF 
energy and the exact energy per site of a fully polar- 
ized state, Ep/N = 2S'^{Ji + J2) (here with S = 1/2). 



Such a point J2 



F-CAF 



(TV) appears, depending on size 
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A^, at a quite precise value of J2/\Ji \ between 0.402 and 
0.4078, as shown in Fig. [TOl The roughly linear depen- 
dence in suggests an extrapolated transition point 
at J2^-^'''' = 0.41|Ji|. 

From the results in this Section we conclude that, on 
the J3 = line, the present SBMFT approach can con- 
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FIG. 7. (color online) Spin correlation function for a system of 
size A'' — lOOx 100 at J2 — 0.75, | Ji| and J3=0 (corresponding 
to the dispersion in Fig. |4]), along both lattice axes. 

netization has a maximum at (±7r, 0). Then Af^(7r, 0), 
shown in Fig. [Sj measures the staggered magnetization 
along the x-direction. In general, M^(7r, 0) and M^{0, tt) 
can be used as order parameters for the CAF phase. 

We now turn to discuss the region 0.4 < J2/IJ1I < 
0.58, under controversy in the literature, where the 
largest sizes studied do not provide a self-consistent so- 
lution. Notice in Figs. [5] and H] that in all converged 
solutions there is no signal of an exotic phase but clear 



J2 > ^2°° = 0.58 I Ji I, this can be tested even in the ther- 
modynamical limit. For lower J2, although systems of 
limited size can be solved, we find CAF observablcs un- 
til the transition to the F phase. A linear extrapolation 
of the CAF ground state energies suggests a direct first 
order transition to the F phase at J2~'^^^ = 0.41 1 Ji| in 
accordance with Ref. 1241. 



IV. EFFECTS OF NEXT TO NEAREST 
NEIGHBOR FRUSTRATION J3 

In this section we analyze the influence of third neigh- 
bors antiferromagnetic couplings J3 on top of the CAF 
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FIG. 9. (color online) Ground state energy in the CAF phase, 
for the same system sizes shown in Fig. [S] Smaller sizes 
(lower panels) show energies obtained for couplings down to 
J2/IJ1I ~ 0.4, in full agreement with extrapolated energies 
from the region J2/IJ1I > 0.58. Exact ferromagnetic phase 
energy (solid line) is included to illustrate the proposed F- 
CAF phase boundary. 
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FIG. 10. Size dependence of J2 



' (A'^), the crossing point 



of ferromagnetic phase energy and the extrapolated CAF en- 
ergy (smallest sizes not shown here). 



phase, by fixmg J2 ~ \Ji\. We recall that, on this Ime, 
the classical phase diagram in Fig. [5] shows coUinear an- 
tiferromagnetic order for < J3 < 0.25| Ji|, a continuous 
transition to coUinear helicoidal (q, 0) (or (0, q)) order 



for 0.25| Ji| < J3 < 0.5|Ji|, with q decreasing from tt to 
Itt, and a discontinuous transition to a helicoidal phase 
{Q,Q) for J3 > 0.5|Ji|, with Q increasing in a narrow 
window, from 0.41957r to ^ (reaching Q = 0.44667r at 
J3 = \Ji\: the largest value of J3 in the present analysis). 

The quantum case was studied by Sindzingre et ali^ 
by ED in systems up to iV = 36 sites, for positive J2 and 
J3, both up to |Ji|. For J2 > 0.75|Ji| and around the 
classical boundary between coUinear helicoidal and heli- 
coidal phases, the authors find signals of an exotic gapped 
phase, stating that is difficult to conclude its precise na- 
ture because of large and irregular finite size effects. In 
particular, on the line J2 = \Ji\: they find a CAF phase 
fo'' ^^3 ^ 0.35| Ji| and a gapped phase for J3 > 0.35| Ji|. 

We have applied the SBMFT to systems of size N = 
20 X 20, 40 X 40, 60 x 60, 80 x 80, 100 x 100. For < J3 < 
0.4| Ji|, we find persistence of the CAF phase: the disper- 
sion relation remains gapless at commensurate momenta 
k ^ (±7r/2, 0) (or (0, ±7r/2)) and the fiLdependent sus- 
ceptibility has a maximum at (±7r, 0) (or (0,±7r)). This 
phase shows a boundary that barely depends on system 
sizes and can be estimated as j^C'^^-C'^f ~ 0.41 |Ji|. 
This amounts to a shift of 0.16|Ji| with respect to the 
classical value. 

For larger J3 the minima of the dispersion relation 
move to incommensurate values oi k. It gets numeri- 
cally difficult on a finite lattice to determine the exis- 
tence of gapless minima. However, though with less pre- 
cision than in the CAF phase, we find for all studied 
sizes that; immediately above J^'' 
coUinear helicoidal phase with k 
vclops. The available values for q are discrete but, as 
shown in Fig. 1111 the dispersion minima position evolves 
in the same range as the classical ones, simply shifted in 
J J, ^ J3 + G.16| Ji|. For each point in the figure, the 0- 
dependent susceptibility shows a maximum at = [q, 0) 
(or (0,(7)), characterizing coUinear helicoidal magnetiza- 
tion order. 

The coUinear helicoidal phase extends up to Js*"^"^ ~ 
0.56| Ji|. Again, such boundary is almost independent of 
the system sizes. For even larger J3 further numerical 
difficulties show up. Indeed, the dispersion seems to get 
gapless at momenta {Q,Q) with Q in the same narrow 
window found in the classical phase diagram for the heli- 
coidal phase. As the thermodynamical values of minima 
position in such a narrow range may mismatch the avail- 
able momenta for a given finite lattice, it is difficult to 
select the minimum of the dispersion relation amongst 
neighboring points. For this same reason, a gap may 
seem to open but could be just a finite size artifact. A 
refinement was done by choosing different system sizes 
in order to allow for different distributions of lattice mo- 
menta. In general, when a clear minimum is found, it 
corresponds to a gapless situation with ordering angle 
9 = {Q, Q) where 0.42l7r < Q < 0.4837r. In such cases, 
the order parameter M^{Q,Q) remains finite signaling 
a helicoidal phase. Thus, the present method shows no 
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V. CONCLUSIONS 



0.42 0.44 0.46 0.48 0.5 
J3/IJ1I 



0.52 0.54 0.56 



FIG. 11. (color online) Evolution of the incommensurate q 
along the line J2 = |Ji|, in the coUinear helicoidal phase for 
different system sizes. For comparison, we show the classical 
pitch angle i?(J3 — 0.16|Ji[) (i.e. plotted with respect to the 
SBMFT phase boundary Jg^^^"^^). 



To summarize, along the line J3=0 and within SBMFT 
we can confirm the coUinear antiferromagnetic phase for 

7^7 F~CAF n /II I 7 I T?^ 7 F-CAF ^ 7 ^ 7 00 

0.58| Ji| the convergence becomes harder, presumably be- 
cause of dominance of ferromagnetic correlations. How- 
ever, at tractable system sizes we find CAF observables 
until the transition to the ferromagnetic phase. A linear 
extrapolation of the CAF phase ground state energies 
suggests a direct first order transition to the ferromag- 
netic phase at J2^~'~^^^ in good agreement with Ref. 

Along the line J2=| Ji| we have found that the bound- 
aries between the coUinear and incommensurate phases 
are strongly shifted, with respect to the classical case, 
to larger values of J3: J^f^^^-^^" « 0.41 |Ji| and 
j^CH-H ^ 0.56|Ji|. We do not find clear evidence of 
spin-gapped phases within the present approximation. 
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